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Abstract 

We show that there is no one-loop enhancement of the rate for a light neutrino to decay 
into a lighter neutrino plus a majoron, contrary to a recent claim. Thus the light neutrinos 
must satisfy the cosmological bound of having masses less than 35 eV in the singlet majoron 
model, or else violate the constraint imposed by galaxy formation. In the latter case, u T could 
have a mass between 40 keV and 1 MeV, while satisfying all other cosmological constraints. 



The singlet majoron model |1] is one of the simplest imaginable extensions of the standard 
model of particle physics. It consists of adding one or more right-handed neutrinos, and a 
complex scalar field a with lepton number —2. If a has a large vacuum expectation value 
v a , one has a natural explanation for the smallness of neutrino masses through the seesaw 
mechanism 0. In addition there is a massless Goldstone boson x, the majoron, which couples 
primarily to neutrinos with a strength suppressed by powers of v a . Although the majoron 
would be impossible to see directly, it could play a role in cosmological and astrophysical 
settings. 

One of the potential uses for majorons is in evading bounds on neutrino masses. For 
example a tau neutrino mass between 35 eV and 1 MeV would seriously conflict with present 
estimates of the energy density of the universe since the standard model decay u T — > is 
too slow to deplete the population of such neutrinos even by the present time. The rate for 
decay into a majoron rather than a photon can be many orders of magnitude faster, thus 
alleviating the conflict. 

In nonminimal majoron models with several scalar fields, it is certainly true that neutrino 
decay rates become unsuppressed But with only one scalar field, it turns out that the 
matrix of couplings of the majoron between the various light neutrinos is proportional to 
the mass matrix, at lowest order in the small parameter e ~ m Ul ^ c /v a [Q. The transition 
v ~> v 'x only occurs at higher order in the small parameter, which typically makes it as 
irrelevant as the electromagnetic decay. It is therefore an interesting claim || || that the 
one-loop correction to the amplitude for v — > v'x is nonzero at leading order in e and is 
diminished only by typical loop suppression factors. This claim has already been used in the 
construction of at least one model @. 

Here we wish to demonstrate the incorrectness of the claim that loop effects enhance the 
neutrino decay rate. We have done this in two different ways. The first method is to compute 
the one-loop effective Lagrangian for the light neutrinos. Let v denote the vector of three 
light neutrino states, and let fi and g denote the 3x3 tree-level mass and majoron-coupling 
matrices respectively. Taking into account one-loop contributions, denoted by <5/x and 5g, 
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the effective Lagrangian becomes 

1 i 
- -v(n + 6n)u- zxPwig + Sg)v, (1) 

We will show that g + Sg is proportional to fi + 5fi at leading order in e, just like the tree-level 
quantities, so that there is no loop enhancement of the decay rate. 

In the effective Lagrangian approach, only one-particle irreducible (1PI) diagrams are 
included. One must then rotate the full matrix of couplings g + 5g to the basis where the 
mass matrix fi + 8[i is diagonal to obtain the physical decay amplitudes. It is also possible 
to calculate the amplitude directly, including the one-particle reducible (1PR) self-energy 
corrections to the tree level vertex, as in reference [fj. An example of such a contribution is 
shown in Fig. 1. We have repeated this computation and found discrepancies with ref. @. 
(We note that it is crucial to put the internal neutrino v<i in the 1PR diagram on the mass 
shell of the external neutrino V\ rather than at zero momentum, in order to get the correct 
result.) In this approach, the 1PR diagrams perform the function of rotating to the one-loop 
mass eigenbasis. We have verified that the 1PR and 1PI diagrams cancel each other, giving 
a vanishing decay amplitude at lowest order in e. Because the effective Lagrangian approach 
is simpler, we will present our results only for that method. 

Before giving the explicit calculations we observe that our result is not a mere accident, 
but can be shown by a general argument. Consider the tree- level Lagrangian for the singlet 
majoron model, 

- hijLiHPuVRj - fijaP R>i P R p R>j + h.c. - A (Jct| 2 - u£/2) , (2) 

where H,L are the Higgs and lepton doublets, P R — (1 + 7s)/2 and v Ri are the isosinglet 
right-handed neutrinos (written in Majorana form). Now imagine doing the renormalizations 
up to any desired order in the symmetric phase, i.e., for negative v 2 a . The mass terms remain 
zero because they are protected by the lepton number and gauge symmetries and only the 
Yukawa couplings get renormalized. The renormalized couplings should have an analytic 
dependence on so that we can compute one- loop effects in the broken phase simply 
by substituting a positive value for v 2 a into the above results. Then the proportionality of 
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resulting renormalized masses and couplings is manifest just as for the tree level parameters, 
to all orders in the loop expansion. While this argument alone might suffice to prove our 
point, we feel obliged to show explicitly how the correct result can be derived in the broken 
phase, where the computation is more involved. 

The first step is to rewrite the Lagrangian when the Higgs and a fields have acquired their 
VEV's. There is a Dirac mass matrix to = hv/\/2 connecting the left- and right-handed 
neutrinos, where v = 246 GeV is the VEV of the Higgs, and a Majorana mass matrix 
M = y/2fv a for the right-handed neutrinos alone. Assuming the eigenvalues of M are much 
larger than the entries of to, the full mass matrix can be partially diagonalized by 



u l\ ( 1 e \ / v 
-e T 1J\N 



e = mM-\ (3) 



so that v and iV respectively have the light and heavy mass matrices /i = — toM 1 to t and 
M, to leading order in e. We also expand 

<? = -j= (Va + P + ix) ■ (4) 

In constrast to the massless majoron x, P has a mass given by m 2 p = 2\v%. The interaction 
Lagrangian is then 

1 , 1 r T ,[eMe T -eM\ , ^fv\ ^ , n 2 2^2 ,r\ 

from which the Feynman rules are easily deduced. The above-emphasized proportionality of 
the vvx coupling matrix to the light mass matrix /i is manifest: 

9 ^ = -ii. (6) 

The next step is to examine the general structure of the one-loop shifts in the mass and 
coupling matrices. In the complete space of the light and heavy neutrinos, (is, N), the masses 
have the form 

1 / n + 5/i Sm 



m " = 2 k T M + 5M'- (7) 
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which necessitates rediagonalizing the light and heavy states using the transformation 

0=(_^ T 5e = 5mM-\ (8) 

Here and below we have kept only the terms which are of leading order in the combined e 
and loop expansion. (For example, keeping SM in the definition of 5e would be consistent 
only if we were doing a two-loop calculation.) After the rotation (||), p + 5p gets a further 
shift —5m T M~ 1 Sm which is equivalent to a two-loop effect and can be ignored. But there is 
a nonnegligible shift in the matrix of couplings, which upon rediagonalization undergoes 

1 ( 9 + 5g g' + Sg'\ 1 / g + 5g — g'5e T — 5eg' T * 



couplings - 2W r + §glT G + 5GJ • 2 V 

(9) 

where * indicates the terms we are not interested in. Thus what we called 5g in eq. ([I]) is 
actually given by the naive one- loop shift plus terms due to heavy-light mixing at one loop. 
We will henceforth use Sg to denote the 'naive' contribution. From this analysis we see that 
it is necessary only to compute the quantities Sg, 8m, and Sp. These are given, respectively, 
by the Feynman diagrams shown in figures 2, 3 and 4. 

First consider Sg, figs. 2a-b. Although there are additional digrams not shown, with v 
rather than N running in the loop, they are of higher order in e. Adding the contributions 
of fig. 2 evaluated at vanishing external momenta, we obtain the matrix equation 

5g = m 2 p v-hl{M,m p )e T , (10) 
io7T z ar — y z y z 

The factor of m 2 p in the numerator here and below always comes from the difference between 
the propagators of x an d p, or from the fact that the XXP coupling is proportional to 
m 2 p . Counting powers of v a (remember that m p ~ M ~ v a ), we see that Sg is of order 
e 2 ~ (m/M) 2 , and is not proportional to the tree level mass matrix p. If this were the end of 
the story, 5g would cause fast decays among the light neutrinos. But we must also consider 
the mass shifts. 
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The Dirac mass shift Sm is given by fig. 3. A straightforward calculation shows that 



Sm = -m 2 p vfl(M, m p )e T . (12) 

Using 5e = 5mM~ 1 from eq. (H) and g' = —eM fv a from eq. (|5|), we find that 

^tot =Sg- g'5e T - 5eg' T = -8g, (13) 

so the net effect of rotating away the induced heavy-light mixing is to reverse the sign of the 
naive shift in the couplings. 

The final singlet Higgs contribution is the direct shift in the Majorana mass matrix of 
the light neutrinos, fig. 4, which we find to be given by 5fi = v a 5g. Therefore the relation 
between the coupling and mass matrices at one loop and leading order in e is 

9 + dg TO T = , (14) 

that is, they are exactly proportional to each other. Thus they are simultaneously diagonal- 
ized and there is no loop-enhancement of decays. 

From the general argument given above, it should be clear that our results do not depend 
on what kind of interactions are involved, so long as they do not explicitly break the global 
lepton number symmetry. Thus the one-loop effects involving Higgs and weak gauge bosons 
must also preserve the proportionality between masses and couplings. We have also carried 
out this calculation. 

The relevant diagrams are shown in figs. 5 and 6 for Sg and Sfi respectively. In contrast 
to the previous case, simple power counting (taking into account necessary insertions of 
neutrino masses) shows that the electroweak contribution to 5m is 0(e 3 ) and so can be 
neglected. The Z-boson interactions may be written as 

i^ro AM^"< (15) 

using the fact that the vector current vanishes for Majorana fields. In evaluating the diagrams 
it proves convenient to work in gauge with £ = M|/M|, for then the contributions from 
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the real and imaginary part of the Higgs field exactly cancel each other, and we are left with 
only the gauge boson diagrams. We obtain 

^ = ~ A 9 ln e ( 3/ ( M > M -) + (M 2 JM 2 Z )I(M, M H )) e T , (16) 
4 cos^ U w v ' 

where / is given in eq. (fll|) ; notice that the relative contributions of the Higgs boson and 
the three polarizations of the Z are evident. We also find that 

Sg TOT = Sg = - — , (17) 

explicitly showing that the weak interactions also preserve the proportionality between 
masses and couplings that forbids transitions between light neutrinos at 0(e 2 ) in the ampli- 
tude. 

(While the contribution of fig. 7 might at first appear to be relevant, it is 0(e 4 ). The 
subgraph containing Z^-x mixing is of order e 2 q fJL) where g M is the majoron momentum. 
Integrating by parts the effective operator v^^vd^x and using the equations of motion gives 
another factor of the light neutrino mass, /i ~ e 2 .) 

Finally, we examine the tree level decay rate more explicitly than appears to have been 
done elsewhere. With a little effort, one can go to next order in e in the diagonalization of the 
neutrino mass matrix, to find that the off-diagonal block of the heavy-light rotation matrix, 
e in eq. (|3]), should be supplemented by the term —mM~ 1 (m T mM~ 1 + ^M~ 1 m T m)M~ 1 . 
This is enough to determine the masses and majoron couplings for the light neutrinos to 
order e 4 , 

fi = fio~Uno,ee T }, ^7 = -— + { — , ee r ) (18) 
2 v a Vv a J 

where /i = —eMe T . To more easily explore the consequences of this result, let us assume 

there is mixing only among the last two generations, with eigenvalues m Vi and Mi for the 

light and heavy masses. In the mass basis, one can show that the off-diagonal coupling in 

eq. (0) is given by 

5-23 = 7j sin 2a [M^ 1 - M 3 _1 ) (m U2 m„ 3 ) 1/2 {m U2 + m U3 )/v a , (19) 
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where a is the rotation angle for diagonalizing the matrix M~ 1 ^ 2 m T mM~ 1 ^ 2 . (Note that 
this vanishes if the heavy neutrinos are degenerate, as expected, since if M is proportional 
to the unit matrix, the Dirac masses can be diagonalized simultaneously with it, in which 
case there is no mixing or decay.) For the purpose of obtaining cosmological upper bounds 
on v a , we conservatively assume that the angle is large and the eigenvalues Mj have a large 
splitting, so that #23 ~ m v( 2m ul 2 / v ai leading to a decay rate of 



1 m V2 m 

— r\j — 



(20) 



for the process u 3 — > v 2 x- 

If m V3 exceeds the cosmological upper bound of 35 eV for stable neutrinos ||, u 3 must 
decay fast enough so that its relativistic decay products have time for their energy to be 
redshifted away. From eq. (EDf) it is clear that we want m U2 to be as large as possible. 
Demanding that the present energy density of v% particles plus majorons from the decaying 
z/ 3 not exceed the closure value, we obtain a bound involving the decay temperature T D , 

m V2 + (ml 2 +p 2 ) +p<35eV; p = (21) 

where p is the redshifted momentum of the majoron. With the aid of eq. ( |20| ) and the 
time-temperature relationship for a universe dominated by nonrelativistic 1/3's, we get 

m V3 \ 1 I 2 



"° < 500 GeV InivJ • (22) 

assuming m U2 has the optimal value (~ 9 eV) for giving the least restrictive bound. 

If m U3 was smaller than 240 keV, this would imply that v a is below the weak scale of 
246 GeV, a somewhat small value from the perspective of the seesaw mechanism for the 
neutrino masses. Nevertheless with v a = 100 GeV, the Dirac neutrino masses corresponding 
to m Ul ~ 1 eV, m U2 ~ 10 eV, m U3 ~ 40 keV are 0.3 MeV, 1 MeV and 60 MeV, respectively, 
comparable to the smallest lepton and quark masses. If m ua saturated the nucleosynthesis 
bound of 500 keV then v a could be as large as 350 GeV and the neutrino Dirac masses 
could span a slightly greater range, up to 0.4 GeV for the third generation. 
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If the majoron has a gravitationally induced mass ||10|| , whose natural value is around 
1 keV, the above results still hold, except for the fact that now the final density of ^2 will 
be four times its thermal value (due to the decay products of x ~ ¥ v \i v i±)i so that its mass 
should be somewhat smaller. The ensuing bound on v a decreases only slightly 

Although considerations of the cosmological energy density leave room for the relevance of 
light neutrino decays in the singlet majoron model, the growth of large scale structure leads 
to a different conclusion. The decay products of a massive tau neutrino could cause a second 
radiation-dominated era, during which density perturbations grow only logarithmically with 
time. The recent measurements of COBE [|TIJ indicate that any additional period of radiation 



domination would conflict with the necessity for primordial perturbations to have grown into 
galaxies by now. Therefore v T should never have matter-dominated the energy density of 
the universe [|12| , hence the decay temperature of v T must not be much less than its mass 



(more precisely, m Uz /T D < 14; see ref. |13[), leading to the constraint 



For any allowed value of m V3 , this is too small a value of v a to give plausible seesaw masses 
for the light neutrinos. 

In summary, we have given both a simple general argument valid to all orders in per- 
turbation theory, and detailed calculations at one loop in the broken phase, to show that 
decay amplitudes between light neutrinos always vanish at order e 2 in the simplest majoron 
model, where e is the ratio between the scales of Dirac and Majorana masses of the neutri- 
nos. Thus there is no enhancement of the amplitude v —>■ v'x due to loop corrections in the 
simplest majoron model. We have also shown that the tree level decays are too slow to be 
cosmologically relevant, but only if the galaxy formation constraint is considered. 
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Figure Captions 



Fig. 1 One-particle reducible contribution to the decay amplitude. 
Fig. 2a,b Singlet Higgs contributions to 8g. 
Fig. 3 Singlet Higgs contribution to 5m. 
Fig. 4 Singlet Higgs contribution to S/i. 

Fig. 5a,b Electroweak contributions to Sg, using all four combinations of v and iV in 5b. 
Fig. 6a,b Electroweak contributions to S/i. 
Fig. 7 A subleading contribution to 5g. 
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